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ABSTRACT 

In this paper we present the results from the analysis of a sample of 28 gamma-ray 
burst (GRB) afterglow spectral energy distributions, spanning the X-ray through to 
near-infrared wavelengths. This is the largest sample of GRB afterglow spectral energy 
distributions thus far studied, providing a strong handle on the optical depth distribu- 
tion of soft X-ray absorption and dust-extinction systems in GRB host galaxies. We 
detect an absorption system within the GRB host galaxy in 79% of the sample, and an 
extinction system in 71% of the sample, and find the Small Magellanic Cloud (SMC) 
extinction law to provide an acceptable fit to the host galaxy extinction profile for 
the majority of cases, consistent with previous findings. The range in the soft X-ray 
absorption to dust-extinction ratio, Nh.x/Av, in GRB host galaxies spans almost two 
orders of magnitude, and the typical ratios are significantly larger than those of the 
Magellanic Clouds or Milky Way. Although dust destruction could be a cause, at least 
in part, for the large Nh,x/Av ratios, the good fit provided by the SMC extinction 
law for the majority of our sample suggests that there is an abundance of small dust 
grains in the GRB environment, which we would expect to have been destroyed if dust 
destruction were responsible for the large Njj,x/ 'Ay ratios. Instead, our analysis sug- 
gests that the distribution of Nh,x /Ay in GRB host galaxies may be mostly intrinsic 
to these galaxies, and this is further substantiated by evidence for a strong negative 
correlation between Nh.x/Av and metallicity for a subsample of GRB hosts with 
known metallicity. Furthermore, we find the Nh,x/A\t ratio and metallicity for this 
subsample of GRBs to be comparable to the relation found in other more metal-rich 
galaxies. 

Key words: gamma-rays: bursts - gamma-ray: observations - galaxies: ISM - dust, 
extinction 



1 INTRODUCTION 

To unravel the properties of gamma-ray burst (GRB) pro- 
genitors and the fundamental conditions required within a 
galaxy to form a GRB, an understanding of the GRB cir- 
cumburst and host galaxy environment is essential. The faint 
optical magnitudes and large distances of GRB hosts limit 
the amount of information obtained from host galaxy ob- 
servations, and thus studying the spectral properties of the 
afterglow provide the most direct and sensitive method of 
probing the surrounding environments of GRBs. 

Both optical and X-ray spectroscopic observations have 



provided detailed information on the metal abundances 
and column densitie s in GRB local environments (e.g. 
Prochaska et al.ll2007l ; ISavaglio et "ail 120031 : ISavaglio fc Fall 
2004l h and broadband analysis of the GRB afterglow spec- 
tral energy distribution (SED) allows the host galaxy dust 
extinction curve to be well modelled, and thus provides 
a mea sure of the host visual extinction. iGalama fc Wiiera 
(2001) combined these two techniques to compare the soft 
X-ray absorption with the visual dust extinction in the local 
environment of a sample of eight GRBs. In their analysis, 
they found that whereas the distribution of equivalent neu- 
tral hydrogen column density within GRB host galaxies was 
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comparable to that observed in Galactic molecular clouds, 
the measured host galaxy visual extinction was 10-100 times 
smaller than expect ed for GRBs embed ded in Galactic- 
like molecular clouds. IStratta et al.l l|2004) expanded on this 
work, and found the gas column density to dust extinction 
ratio to not only be larger than that of the Milky Way, but 
also to be around an order of magnitude larger than that 
of the gas rich Large and Small Mag ellanic Clouds (LMC 
and SMC, respec t ively) . In both iGalama fc Wiiersl (|200lh 
and IStratta et all |200J) the large column density to visual 
extinction ratios measured in GRB local environments was 
taken to be evidence of dust destruction by the GRB, caus- 
ing the visua l extin ction to decrease. In addition to this, 
IStratta et all l|2004l ) also found that the optical to near- 
infrared (NIR) GRB afterglow data showed little evidence 
of the strong 2175 A absorption feature present in the Milky 
Way extinction law, and less promine ntly in the LMC. In - 
stead, they found the mean SMC or the lCalzetti et aD ()l994l ) 
starburst galaxy dust extinction law to provide the best fit 
to their sample of GRB SEDs, both of which have no 2175 A 
absorption feature. T hese results a re sup port ed by the more 
recent work done bv lKann et alj (2006) and IStarling et all 
|2007l ). In each of these cases the mean SMC, LMC and 
Galactic extinction curves were assumed. However, a range 
in the total-to-selective extinction, Rv = Av/E(B — V), 
which relates the reddening to extinction, and 2175 A bump 
strength is obser ved along different lines-of-sight through 
the Mil ky Way (ICardelli et all Il989h and the Magellanic 
Clouds dGordon et al.ll2003l ). and observations of higher red- 
shift supernovae and quasars (z > 5) also suggest differences 
between the dust extinction propert ies in higher redshift 
galaxies and the local universe (e.g. iMaiolino et~ai] |2004 
iTodini fc Ferrarall200lh . However, the degeneracy that ex- 
ists between the best-fit GRB spectral index and the host 
galaxy's total-to-selective extinction means that to accu- 
rately determine the host galaxy's extinction law proper- 
ties, good quality, broadband data are needed, preferentially 
stretching out into the negligibly extinguished far infrared 
(FIR) wavelength bands. 

In the current era of Swift and rapid-response ground- 
based telescopes, prompt arcsecond GRB positions have pro- 
vided a wealth of high quality, early-time X-ray, ultraviolet 
(UV), optical and NIR data. Accurate soft X-ray absorp- 
tion meas urements are now available for the large fraction 
of GRBs rtCampana et~aD l2006bl : iButler fc KocevskJ 120071 : 
iGrupe et al.1 120071 ; ISchadv et all [2007i ; lEvans et al.1 120091 ). 
and well-sampled, high signal-to noise SEDs are providing 
strong constraints on the best-fit extinction law models (e.g. 
iPerlev et al. 1 120081 ). There have now been some examples of 
GRB host galaxies with the 2175 A absorption feature (e.g. 
GRB 070802, Elfasdottir et al. 2009, Kriihler et al. 2008; 
GRB 080607, Prochaska et al. 2009), as well as GRB host 
galaxies with Rv va lues larger than t he mean SMC, LMC 
and MW values fe.g. lPerlev et al.ll2008l ). a possible indicator 
of grey dust, as suggested to be p resent in some GRB host 
galaxies (e.g. ISavaglio et alj|2003l ). However, such analysis 
on the detailed properties of GRB extinction curves are typ- 
ically still only possible for a handful of well-sampled, bright 
GRBs (e.g. GRB050525A, Heng et al. 2008; GRB 061126, 
Perley et al. 2009; GRB 070802, Elfasdottir et al. 2009, 
Kriihle r et al. 2008; GRB 80607, Prochaska et al. 2009). 

In lSchadv et""aH <|2007l ) we used X-ray and UV/optical 



simulta neous observat i ons ta ken with the X-ray Telescope 
(XRT; iBurrows et al. | l2005al) and Ultravi olet and Opti- 
cal Telescope (UVO T; iRoming et~aD l2005h onboard Swift 
(jGehrels et al.l l2004) to analyse the SEDs for a sample of 7 
GRBs. The dust extinction in the GRB host galaxies was 
modelled on the mean SMC, the LMC and the Milky Way 
extinc tion curves using the parameterisations given in IPeil 
1 19921 ). which cover a range in 2175 A bump strengths and 
Rv values. The SMC and LMC extinction curves were found 
to provide the best-fit model for the majorit y of the sam- 
ple, in agreement with previous studies (e.g . IStratta et al.l 
12004 iKann et all 120061 ; IStarling et al.ll2007l i However, we 
also found that, although the gas-to-dust ratio in Swift GRB 
host galaxies were typically larger than those of the Milky 
Way and Magellanic Clouds, the weighted mean was within 
90% confidence of the Magellanic Clouds and Milky Way 
X-ray absorption to optical extinction ratios. 

In this paper we aim to further our previous work, us- 
ing a larger sample of 28 GRBs, and increasing the wave- 
length range of the afterglow SEDs to better constrain the 
absorption and exti nction within the GRB host galaxy. In 
ISchadv etT al. (2007) Swift data alone were used to produce 
the SEDs, and UVOT data with a rest-frame wavelength 
A < 1215 A were not included in the SED fits in order to 
avoid the absorption caused by the Lyman forest being con- 
fused for dust-extinction. In this paper we now model the 
absorption resulting from the Lyman forest such that all rest 
frame UV data redward of the Lyman edge is included in our 
spectral analysis. Furthermore, we also include additional 
ground-based NIR data if available, further increasing the 
spectral range of the SEDs and the degrees of freedom of the 
spectral fits. This provides better sampled SEDs and extends 
the redshift range within our sample, which was previously 
restricted to z < 1.7 to ensure that the SED modelling was 
sufficiently well constrained within the optical wavelength 
range. 

In fJ21 we present the new, extended GRB sample and 
describe the X-ray, UV/optical and NIR data reduction and 
analysis, and in fJS] we describe the models used to fit the 
data. We present the results of our spectral modelling in §3] 
followed by an analysis of the possible selection effects and 
systematic biases that may be present in our work in fJSl A 
discussion on the implications of our findings is presented in 
fj6] and our conclusions are summarised in iJZ] Throughout 
the paper temporal and spectral indices, a and j3, respec- 
tively, are denoted such that F{v, t) oc v~^t~ a , and all errors 
are la unless specified otherwise. 



2 DATA REDUCTION AND ANALYSIS 

The selection criteria for our sample is that the GRB must 
be long (i.e. Tgo> 2 s, where T90 is the time-interval over 
which 90% of the high-energy radiation (> 15 keV) is emit- 
ted), it must have a spectroscopic redshift measurement, 
have been observed by the XRT and UVOT within an hour 
of the prompt emission, have a peak UVOT v-band mag- 
nitude v ^ 19, and be detected by the XRT and in at 
least three UV-IR filters (UVOT and/or ground based). 
The final requirement is needed in order to provide suf- 
ficient constraints for spectral fitting. A total of 28 Swift 
GRBs satisfied our selection criteria up to and including 
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Table 1. Table listing the 28 GRBs in our sample with their rcdshift, Galactic column density and visual extinction in the line-of-sight 
to the GRB, the corresponding SED epoch, the UV, optical, and NIR band passes included in the GRB afterglow SED, and the rest 
frame coverage of the SED. 
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GRB 070529. By requiring that the GRBs in our sample 
have both a spectroscopic redshift measurement and UVOT 
v-band magnitude v 19 we are introducing a bias against 
highly extinguished GRBs, that occur in very dusty regions 
of their host galaxy and/or along a line-of-sight with high 
foreground extinction. A number of previous studies have 
already shown that subsamples of GRBs with spectroscopic 
redshifts are biased against high obscuration (e.g. Fiore et al. 



2007; Fynbo et al. 2009). Further to this, there is also a se- 
lection effect in the redshift distribution that biases against 
certain redshift ranges that have few prominent absorption 
lines in the observer frame optical bandpass, thus making it 
difficult to acquire an accurate spectroscopic redshift mea- 
surement. However, there is currently no evidence to suggest 
that there is a strong redshift dependence on the environ- 
mental conditions of GRB host galaxies, and therefore, for 
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the purposes of this paper, where we are primarily interested 
in studying the dust and metal contents in the environments 
of GRBs, this selection effect in the GRB redshift distribu- 
tion should not degrade our results. Due to the indepen- 
dence between foreground and host galaxy extinction, the 
bias against high foreground extinction should also not have 
any impact on our results on the GRB host extinction and 
absorption properties, and we therefore need to only worry 
about the selective effects introduced by large host galaxy 
extinction, which we explore in detail in section [5~T1 

In order to measure the level of host galaxy dust- 
extinction and absorption in the GRB optical, UV and X-ray 
afterglows, we produced an SED at a single epoch for each 
of the 28 Swift detected GRBs in our sample, where the SED 
epoch was GRB dependent. All Swift data used to produce 
these SEDs were taken from the UK Swift data archive 1 . 
NIR data reported in refereed papers and GCNs were used to 
extend the afterglow SED to longer wavelengths, where pref- 
erence was given to photometry from ref ereed journals. Fur- 
thermore, unlike in ISchadv et all l|2007h . UVOT data with 
rest frame wavelengths A < 1215 A were also used. Absorp- 
tion at these wavelengths caused by the Lyman for est was 
modelled using the work described in iMadaul i|l995l ). which 
provides a statistical estimate of the number of intervening 
absorption systems in the line-of-sight as a function of red- 
shift and column density, and thus opacity of the Lyman 
forest as a function of wavelength. 

The epoch of the SED was chosen to minimise the total 
amount of interpolation required for each UV, optical and 
NIR photometric data point used in the SED. Since XRT 
and UVOT observations are taken simultaneously, this con- 
dition also limited the amount of interpolation required in 
the X-ray band, in which the GRB aftergl ow is typically de- 
tected for longer than in the UVOT (e.g. lOates et af1l2009l ; 
lEvans et al.ll2009T ). A further condition on the selected epoch 
of the afterglow SED was that there could not be any appar- 
ent spectral evolution in either the UVOT or XRT energy 
ranges during the interval used for photometric interpola- 
tion, as is sometimes observed during the early-time stee p 
decay phase o f the X-ray light cu rve (|Nousek et al.| [2006). 
during flares dFalcone et alj 120071 ). or in t he presence of a 
super nova component (e.g. GRB 060218, ICampana et al.l 
l2006h . The 28 GRBs in the sample are listed in Table [TJ 
together with their spectroscopic redshifts, the Galactic hy- 
drogen column density and visual extinction in the line-of- 
sight to the GRB, the epoch of the SED, the UVOT and 
ground-based niters used in the SED, and the rest frame 
wavelength coverage. 

2.1 UVOT and Ground-Based Data 

The UVOT contains three optical and three UV lenticular 
niters, which cover the wavelength range between 1600 A 
and 6000 A, in addition to a clear white filter th at covers the 
wave length range between 1600 A and 8000 A ijPoole et al.l 
2008). The data available to download at the Swift data 
archive are reduced by the science data centre at God- 
dard Space Flight Center, and photometric analysis can 
be carried out immediately on the level 2 products, which 

1 http:/ /www. swift. ac.uk/swift_portal/ 



are already in sky co-ordinates and aspect corrected. In or- 
der to convert UVOT images into spectral files compati- 
ble with the spectral fitting package, xspec, we used the 
tool uvot2pha (vl.l). The response matrices used for the 
UVOT filters were taken from the Swift/VVOT calibration 
files swu** 200^1 120vl04-rsp, where ** is the code for the 
appropriate filter. 

When ground-based optical or NIR data were avail- 
able to use in the SED, spectral files were produced 
for each filter using the appropriate responsivity curves. 
Cousi ns R and I responsivity curves were taken from lBesselll 
(1990), an d the J, H and if-ba nd res ponsivity curves wer e 
taken from lCohen et al I (|l992al lbh and lBessell et all i|l998h . 
For the Sloan Digit al Sky Survey (SDSS) ugriz filters 
(|Fukugita et al.lll996l ). responsivity curves provided in the 
SDSS data release 6 were used 2 . 

To produce an SED at an instantaneous epoch the 
magnitude of the afterglow at the epoch of the SED was 
measured by interpolating or extrapolating the UVOT and 
ground-based filter dependent light curves to the epoch of 
interest. The spectral files were then set to the extrapo- 
lated/interpolated magnitude measured in each correspond- 
ing filter. 

For the UVOT filter light curves, source photometric 
measurements were extracted from the UVOT imaging data 
using the tool uvotmaghist (vl.0) with a circular source 
extraction region that ranged from 3 — 5" radius, depending 
on the brightness of the source. In order to remain com- 
patible with the effective ar ea calibrations, wh ich are based 
on 5" aperture photometry |Poole et al.ll2008T ). an aperture 
correction was applied where necessary. The background was 
taken from a source-free region close to the target with a ra- 
dius of between 10" and 20". The light curves were then 
binned into groups ATu n /T = 0.1, where ATi,i„ is the time 
interval of the bin, and T is the time since the BAT trigger. 

For the ground-based optical and NIR data, filter de- 
pendent light curves were produced using the data from the 
literature. Both for data taken from refereed publications, 
or from GCNs, which are subject to systematic uncertain- 
ties in absolute calibration, the calibration systematic error 
was added in quadrature to the photometric error on each 
measurement. Where no error was provided, either on the 
magnitude or calibration, an error of 0.3 magnitudes was 
assumed. 

When interpolating or extrapolating each filter depen- 
dent light curve to the SED epoch, the same decay index 
was fit over the same time interval for all the filter light 
curves within each GRB (for both UVOT and ground-based 
data). Both the time interval and decay index used were de- 
termined from the combined UVOT and ground-based light 
curve, where all filters were normalised to the UVOT white 
band, if available, and if not, to the v-band, to produce a 
single light curve. The time interval was chosen such that 
it covered the epoch of the SED and could be well-fitted 
by a power-law, and this ranged from AT/Tsed = 0.5 to 
AT/Tsed = 5, where AT is the time interval used, and 
Tsed is the epoch of the SED. The best-fit decay index 
to this time interval was then used to fit each independent 
filter light curve. Having set the spectral files to the extrapo- 



2 http:/ /www. sdss.org/dr6/instruments/imager/ 
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lated/interpolated magnitude measured in each correspond- 
ing UVOT and ground-based filter, a further 10% systematic 
error was added to each ground-based spectral data point 
to account for uncertainties in the responsivity curves. 



2.2 X-Ray Data 

The XRT is well-calibrated in the 0.3-10 keV energy range, 
and has two primary observing modes: window timing 
(WT), which has a 1.7 ms time resolution and 1-dimensional 
imaging, and photon counting (PC), whic h has a 2.5 s 
time r esolution and full imaging capabilities (|Burrows et al.l 
l2005al ). Both modes have spectroscopic capabilities. All data 
were reduced with the xrtpipeline tool (v0.11.6) using the 
most current XRT calibration files, version 20080509. In 
most cases PC mode data were used, with the exceptions 
being GRB 060206 and GRB 061007, which only had WT 
mode data at the epoch of the SEDs (10 ks and 600 s af- 
ter the BAT trigger, respectively). For PC data, both for 
the spectral and temporal analysis, source counts were ex- 
tracted from a circular region centred on the source with a 
radius ranging from 20 to 64 pixels, where 1 XRT pixel is 
2.36". When the source was piled-u p we fitted the sou rce 
PSF profile with XRT's known PSF ifMoretti et al.ll2005l ') to 
determine the radius at which pile-up becomes important, 
and used an annular extraction region to exclude data within 
this radius, which ranged from 4 to 6 pixels. The background 
count rate was estimated from a circular, source-free area in 
the field of view (FOV) with a radius ranging from 42 to 
64 pixels. For WT mode data, the extraction regions used 
for the source and background were slits positioned over the 
source and in a source free region of the FOV, with lengths 
ranging from 20 to 80 pixels, respectively. For both PC and 
WT data, we used xselect (v2.4 3 ) to extract light curves 
and spectral files from the event data in the energy ranges 
0.3-10 keV, which is the band required for compatibility 
with the current calibration files 4 . The spectral files were 
grouped to ^ 20 counts per energy channel, and the light 
curves were binned into time intervals of AT /T = 0.1. Effec- 
tive area files corresponding to the spectral files were created 
using the xrtmkarf tool (vO.5.6), where exposure maps 
were taken into account in order to correct for bad columns. 
Response matrices from version 10 of the XRT calibration 
files were used for both WT and PC mode data. The spec- 
tra were normalised to correspond to the 0.3-10 keV flux of 
the GRB afterglow at the epoch of the SED. This flux was 
determined from the best-fit power-law decay model to the 
afterglow light curve, in the same way as was done for the 
UVOT and ground-based data. 

3 THE MODEL 

The SEDs were fitted within xspec (vl2.4.0 5 ) us i ng th e 
same spectral models as those used in ISchad~ et all (2007), 
with the exception that in this paper absorption due to the 
Lyman forest is also accounted for. 



For each SED we tried both a power-law and a bro- 
ken power-law to fit the afterglow spectral continuum. In 
the latter case the change in spectral slope was fixed to 
A/3 = 0.5 to correspo nd to the chang e in slope caused by 
the cooling frequency (|Sari et al.lll998l ) lying within the ob- 
served frequency range at the epoch of the SED. In both 
the power-law and broken power-law models, two indepen- 
dent dust and gas components were included to correspond 
to the Galactic and the host galaxy photoelectric absorp- 
tion and dust extinction. The Galactic components were 
froze n to the colum n dens ity a nd reddening v a lues t aken 
from Kalberla et al. ( 2005]) and ISchlegel" et al. re- 
spectively, which although uncertain, in particular for lines- 
of-sight with large Galactic reddening, we found to be typ- 
ically an order of magnitude smaller than the errors on 
the measured host galaxy absorption and extinction val- 
ues. Furthermore, the uncertainty in the Galactic values 
becomes negligible when propagated to the additional er- 
ror on the measured rest frame absorption and extinction 
values. The dependence of dust extinction on wavelength in 
the GRB host galaxy was modelled on the SMC, the LMC 
and the Milky Way (MW) empirical extinction laws using 
the XSPEC tool zdust, which is b ased on the extinction co- 
efficients and extinction laws from iPeil (|l992t ) . The total-to- 
selective extinction was taken to be Rv = 2.93, 3.16 and 
3.08 for th e SMC, L MC and Galactic extinction laws, re- 
spectively (|Peilll992T ). From here onwards we shall refer to 
each of the spectral models as the SMC, LMC and MW 
model, where the name corresponds to the extinction law 
used to describe the dust extinction properties in the GRB 
host galaxy. The equivalent neutral hydrogen column den- 
sity in the host galaxy was determined from the soft X-ray 
absorption, where solar abundances were assumed, and is 
denoted throughout this paper as Nh,x- 

There have been a number of examples where in- 
tervening systems have been det ected in the line-o f- 
sight to GR B s (e.g . GRB 050730; IStarling et all <l2005h; 
iD'Elia et all [|2007h . GRB 05092 2 C; iPiranomonte etaL 
(|2008j) . GRB 060418; lElhson et al l (120061); IVreeswiik et al. 
(|2007h . GRB 070802; lEhasdottir et al.1 (|2009l) '). although in 
the majority of cases, the dominant absorption system has 
been reported as originating at the host galaxy. The largest 
reported absorption from an intervening system to date cor- 
responds to an absorption system at z = 2.077 in the line- 
of-sight to GRB 050922 C, which had a column density of 
N H i = 2.0 x 10 20 cm~ 2 (|Piranomonte et ai1l2008l '), which is 
an order of magnitude less than the host galaxy neutral hy- 
drogen column density. A notable exception to this was in 
the case of GRB 060418, which had a strong foreground ab- 
sorber at z = 1.118 with an estimated lower limit on the 
hydrogen column density of Nhi > 1-7 x 10 21 cm -2 and a 
marginally larger extinction than at the host galaxy. Never- 
theless, the percentage of GRBs with reported strong inter- 
vening systems is small (e.g. Prochter et al. 2006, ApJ, 648, 
93; Sudilovsky et al. 2007), and not including their interven- 
ing systems in our SED modelling is, therefore, unlikely to 
affect the overall results of this paper. 

To model the Lyman-series absorption in the 912- 
1215 A rest frame wavelength range, we wrote a local model 



3 http: / /heasarc. nasa.gov/docs / softwarc/lheasoft / ftools / xselect / 

4 http:/ /heasarc. gsfc.nasa.gov/docs/heasarc/caldb/swift/docs/xrt/ 5 http:/ /heasarc. nasa.gov/xanadu/xspec/manual/ 
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Figure 1. The unabsorbed spectral energy distributions for the 28 GRBs in our sample in units of mjy against Hz. In each figure, we 
plot the best-fit absorption and extinction corrected spectral model (solid lines) and data (black data points), as well as the host galaxy 
absorbed and extinguished spectral model (dashed lines) and data (open data points). The extinction curve used in the fit is labelled for 
each SED, as well as the underlying continuum fit to the SED; cither a power-law (pow) or a broken power-law (bknp). 
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Figure 2. Distribution of GRB host galaxy Nh,x for those 
GRBs in the sample with a soft X-ray absorption system de- 
tected at the host galaxy at the 90% confidence level (solid his- 
togram). For those GRBs with no detected host galaxy absorp- 
tion system (90% confidence), the distribution of Njj^x 3cr up- 
per limits is plotted (dotted histogram). The dashed histogram 
is the expected neutral hydrogen column density distribution for 
GRBs that occur within Galactic-like molecular clouds taken from 
iReichart fc Price! J2002f) . 

for XSPEC, and this was included in our fit of the after- 
glow SEDs. T he model used the prescription provided in 
iMadaul (|l995t ) to estimate the effective optical depth from 
the Lyman-series as a function of wavelength and redshift. 
As well as estimati ng the hydroge n absorption caused by in- 
tervening systems, IMadaul l|l995l ) also determined the error 
on this due to statistical fluctuations in the number of ab- 
sorbing clouds along the line-of-sight. This error was added 
in quadrature to the photometric uncertainty of any optical 
data at rest frame wavelengths blueward of Lya. 

4 RESULTS 

4.1 Distribution of Host Av and Nh,x 

The results from our spectral analysis are summarised in Ta- 
ble 4, and the absorption corrected SEDs are shown in Fig. \T\ 
in units of mjy against Hz. Each SED in Fig[T]shows the host 
galaxy absorbed and extinguished spectral model and data 
(dashed lines and open data points, respectively), as well 
as the best-fit absorption and extinction corrected spectral 
model and data (solid lines and black data points, respec- 
tively). We find that 79% of the sample has a host galaxy 
absorption system detected with 90% confidence, with an 
equivalent neutral hydrogen column density, Nh,x, ranging 
from N H ,x = 8.2 x 10 20 cm -2 to N H ,x = 1.4 x 10 22 cm -2 , 
and a dust extinction system local to the GRB is detected in 
71% of the sample with 90% confidence, with a range in ex- 



Figure 3. Distribution of GRB host galaxy visual extinction for 
those GRBs in the sample with a dust extinction system detected 
at the host galaxy at the 90% confidence level (solid histogram). 
For those GRBs with no detected host galaxy extinction system 
(90% confidence), the distribution of av 3<r upper limits is plotted 
(dotted histogram). The dashed line corresponds to the expected 
GRB host galaxy Ay distribution when selection effects are taken 
into account (see section (ETTJl . 



tinction of 0.03 < Ay < 0.75. We note here that the best-fit 
visual extinction is dependent on the shape of the e xtinction 
law used to fit the data, and lCardelli et ail l|l989l ) found a 
linear correlation between the total-to-selective extinction, 
Rv , and the amount of U V extinction along different lines- 
of-sight within the Milky Way. For those GRBs best-fit by 
a broken power-law continuum, the additional free param- 
eter in the fit introduces some degeneracy between the lo- 
cation of the spectral break and the total-to-selective ex- 
tinction, Rv, and in these cases, extinction laws with larger 
total-to-selective extinction, Rv, than those considered in 
our analysis could result in larger best -fit extinction values 
than those listed in T able 4 (e.g. see IWatson et alj 120061 ; 
lElfasdottir et al1l2009l ). If we, therefore, only consider those 
results from the GRBs best-fit by a power-law continuum, 
we find that the largest host galaxy visual extinction mea- 
sured in the sample is still Av = 0.75. 

Our measured distribution of host galaxy logarithmic 
column densities for GRBs with a host galaxy x-ray absorp- 
tion system detected at 90% confidence is shown in Fig. [2] 
in units of cm -2 (solid histogram), and has a mean of 21.7 
and a standard deviation of 0.3. For those GRBs that do not 
have a host absorption system detected at 90% confidence, 
the Nh,x 3(7 upper limit distribution is shown by the dot- 
ted histogram. As well as plotting the measured Nh,x dis- 
tribution, we also show the expected Nhi column density 
distribution in the line-of-sight to GRBs within Galactic- 
like m olecular clouds (dashed histogram; IReichart fc Price! 
2002). The expected Nhi distribution will only equate to 
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an Nh.x distribution in the case of solar metallicity GRB 
host galaxies. However, the host galaxies of GRBs typically 
have s ub-solar metallicities that range from 1/100 *'' solar to 



solar ( Prochaska et al. 20071 ; ISavaglio et al. 20091). C o mbin - 
ing the metallicity estim ates from iProchaska et all (|2007l ) 
and ISavaglio et all (2009) gives a distribution with median 
value of l/4 th solar, covering over two orders of magnitude. 
The effect of this on the expected GRB host galaxy Nh.x 
distribution shown in Fig. [2] would be to broaden it to the 
left of the figure, down to values of log Nh,x= 19, and shift 
the logarithm of the expected peak Nh.x value down to 
21.2. However, due to the small number of GRB host galax- 
ies with accurate metallicity measurments, the GRB host 
metallicity distribution is poorly known. We therefore feel 
that for the time being it is not feasible to apply a correction 
to the expected Nhi distribution shown in Fig.[2]to convert 
it into an accurate expected Nh,x distribution. 

Despite the above caveat with the expected Nh,x dis- 
tribution, it appears from Fig.[2]that those GRBs with host 
galaxy column densities at the very low and very high end 
of the distribution are missing from our sample. The lack 
of GRBs in our sample with measurements of host galaxy 
Nh.x smaller than ~ 10 21 cm" 2 is likely to be due to the 
sensitivity limit of the XRT at measuring soft X-ray absorp- 
tion. In which case those GRBs with expected host galaxy 
Nh.x smaller than ~ 10 21 cm -2 may be accounted for by 
those GRBs in our sample with just upper limits on their 
host galaxy soft x-ray absorption. The absence of GRBs in 
our sample with host galaxy Nh.x J> 10 22 cm~ 2 , on the 
other hand, is likely to be a consequence of our selection ef- 
fects. By selecting only those GRBs with observed v < 19 we 
exclude highly extinguished GRBs from our sample, which 
would simultaneously exclude those with large X-ray absorp- 
tion values. We shall address the effect that this bias has on 
our results later, in section[5] Nevertheless, there is generally 
a good agreement between the two histograms, suggesting 
that our selection effects on the distribution of host galaxy 
column densities in the line-of-sight to GRBs is not highly 
significant. 

The measured optical extinction distribution for GRBs 
with a host dust extinction system detected at 90% confi- 
dence is shown in Fig. [3] (solid histogram), and has a mean 
Av of 0.3 with a standard deviation of 0.2. The dotted his- 
togram shows the Av 3<r upper limit distribution for those 
GRBs with no host extinction system detected at 90% con- 
fidence. In this figure we also show the expected GRB host 
galaxy Av distribution when selection effects are taken into 
account (dashed histogram), the details of which we describe 
in section [5Tl 



4.2 Host Galaxy Extinction Curves 

The extinction properties of dust are dependent on the 
dust composition and grain size distribution. The extinc- 
tion curve models that best-fit the data, therefore, provide 
information on the dust properties of the GRB circumburst 
medium. For 18% of the sample no distinction could be made 
in the quality of the fits between extinction curve models, 
and this is a consequence of the lack of absorbing dust in the 
local environments of these GRBs, as well as the low signal 
to noise of some of the data. Evidence of this is provided 
in Fig. [4] where we plot the GRB host galaxy Nh.x against 



j4^in log — log space, using the best-fit values from the SMC 
(top panel), the LMC (middle panel), and the MW spectral 
model (bottom panel). Those GRBs where no distinction 
was possible between dust model fits are shown as grey cir- 
cles, all of which only have upper limit measurements on the 
host galaxy extinction value. 

For those GRB SEDs best-fit by a broken power-law, a 
degeneracy exists between the total-to-selective extinction, 
Rv, and the location of the spectral break, thus limiting 
our knowledge of the shape of the host galaxy extinction 
law in these cases. Where the optical and X-ray afterglow 
emission lie on the same power-law component, however, 
the fewer number of variable parameters provides a greater 
handle on the shape of the host galaxy extinction law. Of 
those GRBs for which a distinction between model fits was 
possible, 21 were best-fit by a power-law spectral model, the 
large majority of which have a rest-frame wavelength range 
that safely covers the location of the 2175A absorption fea- 
ture, and all have data blueward of 1500A in the rest-frame. 
The data available for this subset of GRBs therefore covers 
the wavelength range over which the three extinction laws 
modelled in this paper can be most effectively discriminated 
between, and we therefore use this subset of GRBs to study 
the distribution in best-fit extinction laws. 

We find that the SMC extinction curve provided the 
best-fit in 56% of cases, and the LMC and MW extinction 
curve both provided the best-fit to 22% of cases. Of the four 
GRBs where the MW extinction law provided the best fit to 
the SED, three have optical data that span the wavelength 
range of the 2175A absorption bump at the rest frame of the 
host galaxy (GRB 050802, GRB 050922C, GRB 070110), 
making the detection of the 2175 A feature possible. Never- 
theless, in all three cases the SMC and LMC spectral models 
also provide acceptable fits, and we, therefore, cannot claim 
a robust detection of the 2175 A absorption feature in the 
host galaxy of these three GRBs. The question of how typical 
the Milky Way absorption feature is in GRB host galaxies 
is beyond the scope of this paper, and an issue that we are 
looking to investigate in future work. That the host galaxy 
dust extinction properties for the majority of our sample are 
best-fit by the SMC extinction law is consistent with sev- 
eral previous stu d ies in this field (e.g . IStratta et al. 1 12004 
iKann et al.1 120061 ; IStarling et all I2007T ). It is not yet clear 
what is responsible for the Milky Way absorption feature at 
2175 A, although small carbon acious dust grains ar e thought 
to play an important role (e.g. iDraine fc LeeHl984T ). suggest- 
ing that such grains are not typical in the environments of 
GRBs, at least not once the GRB has occurred. 



4.3 Nh.x/Av Ratio In GRB Host Galaxies 



In Fig. |4j the dashed lines in each panel represent 
the mean hydrogen column densi t y to extinction ratios 
in the SMC (t op; iMartin et all ll989l ). LMC (middle; 



KoornneeJI 19821; iFitzpatric k 1985) and Milky Way (bottom; 



Predehl fe Schrmttlll995l )T which have been converted from 
Nhi/Av to an Nh,x/Av ratio relating to the column den- 
sity that would be measured from X-ray observations of each 
of these galaxies if solar abundances were assumed. We did 
this by assuming a me tallicity of 0.25 sol ar for the SMC, and 
0.5 solar for the LMC l|Wood et al.ll 19981 ). and the N H ,x /A v 
ratio was then a fraction 0.25 and 0.5 the Nhi/Av ratio 
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Figure 4. Host galaxy Ay against Nh.x in log — log space, taken 
from the spectral results from the SMC model (top panel), the 
LMC model (middle panel) and the MW model (bottom panel) 
spectral fits. The dashed and dotted curves are the Nh.x / Ay 
ratios and lcr deviations, respectively, for each corresponding en- 
vironment, where we converted Nh to an X-ray equivalent Nh,x 
value assuming a metallicity 0.25 and 0.5 solar for the SMC and 
LMC, respectively (see text for details). The open squares corre- 
spond to those GRBs that were also used in the ISchadv et al.l 
(2007) sample, and the grey filled circles correspond to those 
GRBs for which no distinction can be made between the goodness 
of fit of the three spectral dust models. 

for the SMC and LMC, respectively. As well as the mean, 
we also show the root-mean square deviation (dotted lines) 
of the sample used to derive the mean SMC, LMC and 
Milky Way Nhi/Av ratios, also converted into the equiv- 
alent Nh.x /Ay ratio. The sub sample of GRBs that were 
analysed in lSchadv et al.l (|2007| ) are shown as open squares. 

From Fig. Uthere does not appear to be any strong cor- 
relation between the dust and gas column density in GRB 
host galaxies. Using only those GRBs with an extinction and 
absorption system detected with 90% confidence, a spear- 
man rank test between the best-fit Av and Nh,x measure- 
ments from the SMC, the LMC and the MW models indi- 
cates a weak correlation at the la level, with coefficients 
0.39, 0.37 and 0.48, respectively. It is also notable that most 
of the data points lie to the right of the dashed lines, corre- 
sponding to Nh,x/Av ratio s that are larger tha n those of 
the SMC, LMC and MW. In lSchadv etldli|2007h the GRB 
Nh.x I Av ratios were also typically larger than those of the 
MW and Magellanic Clouds, although they were still con- 
sistent at the 68% confidence level with the SMC, LMC and 
MW Nh.x I Ay ratios. 

The larger sample used in this paper shows a spread 
in gas-to-dust ratios that covers nearly two orders of mag- 
nitude. This is better illustrated in Fig. [5] where we plot 
the Nh,x/Av ratio from the SMC (top), LMC (middle) 
and MW (bottom) model fits for each GRB in the sam- 
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Figure 5. Ratio of rest frame Nh.x to Ay resulting from the 
SMC (top panel), LMC (middle panel) and MW (bottom panel) 
extinction spectral models, where the results from the best-fit 
model between either the power-law or broken power-law fit are 
shown. The dotted lines correspond to the SMC, LMC and MW 
empirical Nh.x / Ay ratios, and the dashed lines indicate the 
mean Nh.x I Ay value from each set of extinction models fitted. 

pie with a host galaxy absorption and/or extinction sys- 
tem detected with 90% confidence. In each case the data 
points shown are the results from the best-fit model to the 
continuum, i.e. either a power-law or broken power-law fit 
to the SED. The dashed lines are the mean Nh.x I 'Ay ra- 
tio of our sample determined from our spectral analysis for 
each of the extinction curve models fitted. These correspond 
to {Nh.x I Ay) = 3.3 x 10 22 cm -2 , 3.4 x 10 22 cm -2 and 
2.1 x 10 22 cm -2 , which is a factor of 8.3, 9.7 and 11.7 larger 
than the mean Nh.x I Ay ratios measured in the SMC, LMC 
and MW (dotted lines) in the top, middle and bottom panels 
in Fig. [5] respectively. The standard deviation of the data 
about the dashed lines is 2.8 x 10 22 cm -2 , 5.9 x 10 22 cm~ 2 
and 1.8 x 10 22 cm" 2 for the SMC, the LMC and the MW 
spectral model results, respectively. 

5 SELECTION AND SYSTEMATIC EFFECTS 
5.1 Selection effects in Ay 

Those GRBs located in very dusty regions are less likely to 
be detected at optical wavelengths than those GRBs with a 
small host galaxy extinction, and could, therefore, be miss- 
ing in our sample, as was pointed out in section |4~T1 A pos- 
sible indication of this selection effect is the small number 
of GRBs that have Ay> 1 in all three panels of Fig. [4] 
To ascertain better the impact of these selection effects on 
our results we compared our measured distribution of host 
galaxy Ay and Nh.x with the distribution resulting from a 
Monte Carlo simulation of 1000 GRBs with host galaxy neu- 
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Figure 6. Our sample distribution of the extinction corrected 
GRB v-band magnitude. 



Table 2. Results from two-dimensional KS test between the ob- 
served host galaxy ATjj x an d Ay distribution and the distribu- 
tion simulated in eight separate Monte Carlo simulations. The 
simulated host galaxy Ay values are calculated from a randomly 
selected host galaxy neutral hydrogen column density, Njj x> as- 
suming an Nfl x I Ay ratio selected by random from a Gaus- 
sian distribution. The Gaussian Nh,x / Ay distribution assumed 
in each Monte Carlo simulation executed have different pairs 
of mean Nh,x / Ay ratio and Nh,x / Ay standard deviation, a, 
given in columns 2 and 3 respectively. Columns 4 and 5 list the 
KS statistic, D, and KS probability, and in column 7 we give the 
fraction of GRBs rejected by our selection criteria due to having 
an observed magnitude v > 19. 



Model 


mean 


(j 


KS stat. 


KS 


Fraction of 


number 


N H ,x/A v 




D 


Prob 


sample with 




(10 22 ) 


(10 22 ) 






v > 19 


1* 


0.4 


1.0 


0.622 


0.000 


57.20% 


2 


0.4 


2.0 


0.430 


0.007 


57.30% 


3 


1.3 


1.0 


0.311 


0.095 


31.10% 


4 


1.3 


2.0 


0.306 


0.108 


35.20% 


5 


4.0 


1.0 


0.344 


0.047 


85.80% 


6 


4.0 


2.0 


0.386 


0.188 


19.80% 


7 


13 


1.0 


0.581 


0.000 


7.90% 


8 


13 


2.0 


0.587 


0.000 


11.30% 



SMC mean Nh.x / Ay and standard deviation. 



tral hydrogen column densities, Nh,x, and host visual ex- 
tinctions, Ay taken at random, where a Gaussian Nh.x / Ay 
distribution was assumed. 

We selected at random an equivalent neutral hydrogen 
column density, Nh,x, from the expected GRB host galaxy 
Nh,x distribution shown in Fig. [2] For a given Nh,x, the 
host galaxy visual extinction, Ay, was then determined by 
assuming a random column density to visual extinction ra- 
tio, Nh,x I 'Ay , taken from a Gaussian distribution with a 
mean Nn,x/Ay ratio and standard deviation equal to that 
of the SMC. We chose to use the SMC N H ,x/A v distribu- 
tion since the majority of our sample are best-fit by an SMC 
host galaxy extinction law. In addition to selecting a host 
galaxy Nh,x and Ay, we also selected at random a redshift 
with z ^ 4 from the known Swift distribution, and a Galac- 
tic visual extinction and extinction corrected GRB appar- 
ent v-band magnitude from our sample distribution, shown 
in Fig. [U With these parameters we could then determine 
the extinguished v-band magnitude that would be observed 
in each case, thus allowing us to calculate the fraction of 
simulated GRBs with observed v-band magnitudes v > 19 
that would thus be rejected by our sample selection crite- 
ria. Using only those generated GRBs with observed v-band 
magnitudes v < 19, we then performed a two-dimensional 
Kolmogorov-Smirnov (KS) test between the simulated and 
measured Nh,x and Ay data sets, and found that they 
had less than a 1.2 x 10~ 5 probability of coming from the 
same parent population. The generated distribution of Ay 
against Nh.x is shown in log-log space in the top panel of 
Fig. [7] (small, open circles and triangles), together with our 
measured Ay, Nh,x values (solid circles). The poor agree- 
ment between the generated and observed data points can 
be clearly seen in this plot, with the generated data points 
typically lying above the observed data points. 

We, therefore, re-ran our Monte Carlo simulation, but 



this time we used a Gaussian Nn,x/Ay distribution with 
a standard deviation two times that observed in the SMC, 
as well as trying a distribution with a mean Nh,x I Ay ratio 
0.5 dex, 1.0 dex and 1.5 dex larger than that of the SMC, 
each with a standard deviation equal to and two times that 
observed in the SMC. The mean Nn,x/Ay ratios and stan- 
dard deviations of the Gaussian distributions used in each 
of our Monte Carlo simulations are listed in Table [2] along 
with the results from a KS test between our measured Nh,x 
and Ay distribution and the simulated Nh,x and Ay data. 
In the last column of Table [2] we also give the fraction of 
data points from each of the Monte Carlo simulations re- 
jected by our selection criteria as a consequence of having 
an extinguished v-band magnitude v > 19. 

Of all the Nh.x /Ay Gaussian distributions that we 
tried, we found that the two distributions of Nh,x / Ay with 
mean 0.5 dex larger than that of the SMC (models 3 and 4), 
and the two with mean Nh.x I Ay ratio 1.0 dex larger (mod- 
els 5 and 6) all produced samples of Nh.x and Ay values 
that were consistant with our samples of host galaxy Nh,x, 
Ay measurements. In each of these cases the KS probability 
of the simulated and the measured data sets coming from 
the same parent population was at least 4%. The simulated 
data sets that used an Nh,x/Av distribution with the same 
mean Nn,x/Ay as that of the SMC (models 1 and 2), and 
a mean Nn,x/Ay ratio 1.5 dex larger than that of the SMC 
(models 7 and 8), had a KS probability of less than 1% of 
being consistent with our observed sample. 

Another result of our Monte Carlo simulations is that 
they also provide the fraction of simulated GRBs rejected 
by our selection criteria. In a detailed analysis on 14 'dark 6 ' 
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Figure 7. Best-fit host galaxy Ay against Nh X ln log-log space 
for the sample of GRBs analysed in this paper (large circles), plot- 
ted together with a sample of randomly generated host galaxy Ay 
and Nh.x values. In the top panel a Gaussian Nh xl Ay distri- 
bution with the same mean and standard deviation as observed 
in the SMC was assumed for the simulated data, and in the bot- 
tom panel a Gaussian distribution was again assumed, but this 
time with a mean Njj,x I Ay ratio 0.5 dex larger than that of 
the SMC, and two times the standard deviation. In both panels, 
those generated data with corresponding observed GRB v-band 
magnitudes v < 19 are plotted as stars, and those with v > 19 
are shown as small open circles. The mean SMC Nh,x I Ay ra- 
tio (dashed line) and lcr deviation (dotted line) are shown as a 
reference. 



GRBs, iPerlev et all (|2009h estimated that at least 45% of 
those GRBs with an apparent peak magnitude v > 19 were 
'dark' as a result of dust extinction. They also estimated that 
at least 20% of all Swift GRBs have a host galaxy Av > 0.8, 
and at least 10% have Ay > 2.5. The host galaxy Ay dis- 
tribution produced b y our model 4 is in closest agreement 
with the results from lPerlev et ail l|2009l ). generating a host 
galaxy visual extinction Ay > 0.8 29% of the time, and 
A v > 2.5 12% of the time. Of all our Monte Carlo sim- 
ulations, the fraction of GRBs rejected in model 4 by our 
selection criteria (35%) was also the most consistent with 
the fraction of GR Bs with observed v > 19 estimated by 
IPerlev et al.l ([2009) . Given the similarity between the host 
galaxy Ay and observed v-band di stribution produce d by 
our model 4 and the estimates from lPerlev et al.l (|2009f ). we 
take the results from model 4 to be the most representative 
of the true host galaxy Nh.x and Ay distribution, and thus 
use these to quantify our selection effects. 

In Fig. [7J we plot the distribution of simulated host 



6 Dark here refers to GRBs with no detected optical afterglow or 
with an optical flux that is significantly dimmer than expected 
from the observed X-ray afterglow. 



galaxy Nh.x and Ay values resulting from model 4 in the 
bottom panel, together with model 1 in the top panel for 
comparison. Those GRBs with v < 19 are plotted as open 
stars, and those with v > 19 are shown as small open circles. 
Also shown are the best-fit host galaxy Nh,x and Ay values 
for our observed sample of GRBs (large circles), as well as 
the mean SMC Nh.x/ Ay ratio (dashed line) and standard 
deviation (dotted line). The distribution of simulated data 
points for both GRBs with v > 19 and v < 19 shown in 
the bottom panel of Fig. [7J suggest that even when selection 
effects are taken into account, the majority of GRBs con- 
tinue to have host galaxy Nh.x/ Ay ratios that are larger 
than the SMC distribution. In fact, ~ 80% of simulated 
GRBs have Nn.x/Ay host galaxy ratios larger than the 
mean SMC value, and cover nearly four orders of magni- 
tude in Nh.x/Av. We therefore conclude that the results 
from our analysis that the distribution in the Nn.x/Ay ra- 
tio in GRB host galaxies is broad and typically larger than 
those of the Milky Way and Magellanic Clouds applies even 
when dust selection effects are taken into account. 



5.2 Systematic effects in measuring Nh,x 

Any spectral curvature present within the X-ray energy 
range, such as when the cooling frequency, v c , or the prompt 
emission peak energy, E p k, lies within the X-ray band, can 
result in an over-estimation of the measured Nh,x if it is not 
taken into account in the spectral model fit. By fitting all our 
SEDs with a power-law continuum as well as with a spec- 
tral break corresponding to the cooling frequency, we tested 
for the possibility of v c lying within the observing band and 
used the results from the best-fit model. This, therefore, re- 
moves the probability that any spectral curvature resulting 
from v c lying within the X-ray band was incorrectly inter- 
preted as soft X-ray absorption at the host galaxy. However, 
there may still be spectral curvature within the X-ray band if 
there was ongoing X-ray emission from the GRB during the 
time interval over which t he spectrum was ex t racted . From 
the analysis of 59 GRBs, iButler fc Kocevskil ([2007) found 
such a prompt emission contribution out to a maximum of 
Tgo+10 4 s, after which none of the GRBs in their sample 
showed evidence of spectral evolution. The systematic effect 
that such spectral curvature would have on our SED analy- 
sis, therefore, only applies for those GRBs in our sample for 
which we produced SEDs at an epoch earlier than T+ 10 4 s. 
For the twelve GRBs in our sample f or which this appli es, we 
fit the hardness ratio for each GRB (|Evans et alj |2009) from 
the start of the time interval over which the XRT spectrum 
was extracted onwards, and found no evidence for spectral 
evolution over the time interval fitted. We therefore do not 
believe that spectral curvature in the X-ray band is system- 
atically overestimating the host galaxy Nh.x, in agreement 
with the results found by iNardini et all ([2009), who con- 
cluded that intrinsic curvature in the spectrum could not 
be considered as a general solution for the large GRB host 
galaxy N H ,x- 
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Figure 8. Host galaxy Av resulting from spectral modelling of 
the SED where only those optical/NIR data redward of 1215 A 
in the rest frame were fitted, against the best-fit Ay from fits to 
the complete SED. The spectral results for an SMC, LMC and 
MW host galaxy extinction law are plotted in the top, middle and 
bottom panels, respectively, and the dashed lines correspond to 
where there is no difference between the visual extinctions plotted 
along the x and y-axes. 



5.3 Systematic effects in measuring Av 

5.3.1 Hydrogen absorption versus dust attenuation 

If our Lyman forest model is overestimating the attenuation 
of UV light, this would result in a systematic underestima- 
tion of the extinction of UV light from dust. Therefore, to 
test the accuracy with which we model the Lyman-series ab- 
sorption in our spectral analysis, we re-fitted the SEDs using 
only those optical/NIR data redward of 1215 A in the rest 
frame, where Lyman-series absorption no longer applies. In 
Fig. [8] we plot the best- fit dust-extinction yielded when op- 
tical data blueward of 1215 A were excluded, against the 
best-fit host galaxy dust-extinction values determined from 
the method outlined in section [3] The data points plotted 
in the three panels of Fig. [8] correspond to the best-fit dust- 
extinction values resulting from the SMC (top) , LMC (mid- 
dle) and MW (bottom) models, where the dashed lines cor- 
respond to where there is no difference between the visual 
extinctions plotted along the x and y-axes. In all three pan- 
els the data points are evenly distributed about the dashed 
line. There is, therefore, no evidence to suggest that our 
modelling of the optical depth from the Lyman-series is re- 
sulting in a systematic effect on our best-fit Av values. The 
larger scatter around the dashed line in the bottom panel of 
Fig. [8] is due to the typically poorer fits given by the MW 
model compared to the SMC and LMC spectral models. It 
is worth pointing out that Fig. [8] also suggests that the lack 
of a host galaxy neutral hydrogen absorption component in 
our SED model is not overly affecting our results. 



5.3.2 Dust Extinction Cross-Section 

By only considering the mean SMC, LMC and MW ex- 
tinction laws in our spectral analysis, we may be introduc- 
ing another systematic effect on our results. The amount 
of UV, optical and NIR radiation that is extinguished by 
dust is dependent on th e density of dust in the environment 
l|Perna fc Lazzatl 120021 ). as well as the grain size distribu- 
tion, the grain morphology, and the chemical composition, 
all of which influence the d ust extinc tion cross-section as a 
function of wavelength fe.g |Peilll992l ). The measured visual 
extinction, Av, will therefore d epend strongly on t he ex- 
tinction law fit to the data (e.g. ICardelli et al.| [l989l. This 
is illustrated by the differences in the best-fit Av values 
measured between the spectral fits to the GRB SEDs. The 
total-to-selective extinction, Rv, in the local environment of 
the GRB may be lar ger than the three mean values that we 
have considered (e.g. lPerlev et al]|2008h . representative of an 
environment with a dust size distribution skewed to larger 
grains. This would produce an extinction curve that is flat- 
ter in the NIR wavelength range than the SMC, LMC and 
MW extinction laws. Such an extinction law could account 
for the discrepancy that exists between the small amount o f 
reddening observed in GRB SEDs ("e.g. lStratta et al.ll2005r i. 
and the larger host galaxy dus t column densities d erived 
from dust depletion studies (e.g. ISavaglio et al1 l2003). Such 
an extinction law is caused by a dust distribution skewed 
towards larger grains, which may result if there is ongoing 
dust destruction, or if the GRB itself preferentially destroys 
the smaller dust grains during its initial outburst. In the 
case where the dust in the GRB surrounding environment 
has such a grain size distribution, modelling the optical af- 
terglow SED with an SMC extinction law would underes- 
timate Av, since the steepness of the SMC extinction law 
over the NIR, optical and UV range, would yield a smaller 
value of Av for the same amount of UV extinction. 

An estimate of the host galaxy visual extinction that 
does not require knowledge of the host galaxy extinction 
law is provided from GRB optical spectral analysis, where 
measured metal column densities and dust depletion mod- 
els are used to estimate the fraction of metals locked up 
in grains. From the ana lysis of three GRB optical spec- 
tra, ISavaglio et al.l (|2003l ) measured a mean host galaxy vi- 
sual extinction of (Av) ~ 1.0, which is several times larger 
that th e typical host galaxy extinction values m e asured from 
SEDs l|Stratta et alj|2004 Istarling et al.ll2007l ; iKann et al.l 
2006), such as in this paper. One explanation for this dif- 
ference in the visual extinction estimates between SED and 
optical spectroscopic analysis could be the presence of grey 
dust, which produces a flat extinction law. In such a case 
the near-uniform dust extinction across the optical and UV 
wavelength range would leave the observed optical/UV spec- 
tral slope relatively unchanged, and this may thus result 
in an underestimation of the best-fit Av when fitting the 
SED. However, in the case where the GRB SED is best- 
fit by a single power-law component, the reduced effect of 
dust-extinction on the X-ray and NIR bands allows the un- 
derlying spectral index to be well pinned, and in such a case 
grey extinction should be well detected. Of the 28 GRBs in 
our sample, 21 were statistically better-fit by a power-law 
spectral model, and the majority of our sample should there- 
fore have well-determined measurements of the host galaxy 
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extinction. Furthermore, the results from our MW extinc- 
tion model, which has a natter extinction law than the SMC 
and LMC, and hence yields a larger Av value for a given 
amount of dust absorption in the UV, still gives Nh,x/Av 
ratios that are around an order of magnitude larger than 
that observed in the MW. 

Another important point to consider is that extinction 
estimates from optical spectroscopic analysis require cer- 
tain assumptions to be made on the local environment of 
the GRB, such as the dust depletio n pattern and ionisa- 
tion state of the surrounding gas. In ISavaglio et ail |2003), 
both the GRB host galaxy dust depletion chemistry and the 
Nhi to Ay ratio were assumed to be the same as in the 
Milky Way. In particular, if an LMC or SMC Nm/Ay ra - 
tio were assumed, the Ay estimates in ISavadio et all (|2003l ) 
would decrease by a fac tor of more than 3 and 8, respectively. 
IProchaska et alj (|2007f ) estimated the GRB host galaxy ex- 
tinction for a sample of GRBs using the metal column den- 
sities that they measured in the GRB optical spectra, but 
they assumed an SMC gas-to-dust ratio, and for those GRBs 
where a visual extinction estimate was possible, they esti- 
mated a maximum extinction value of Ay= 0.18. 

Finally, even if we w ere to adopt a n ext inction of 
Av ~ 1.0, as estimated by ISavaglio et alj ([2003), it is still 
considerably smaller than would be expected in environ- 
ments with similar Nh,x/Av ratios as in the Milky Way 
and Magellanic Clouds given the mean GRB host galaxy 
column density, Nh,x- GRB host galaxy extinctions would 
need to be around an order of magnitude larger than the 
values that we measure in our analysis in order to be consis- 
tent with the Nh,x/Av ratios measured in the Milky Way 
and Magellanic Clouds (see Fig. [5]). 



6 DISCUSSION 



In ISchadv etafl (|2007f) we proposed that the general con- 
sistency between GRB host galaxy gas-to-dust ratios and 
those of the Milky Way and Magellanic Clouds was evidence 
that the level of photo-ionisation caused by the GRB was 
of the same order as the amount of dust destroyed by th e 
prompt emission, as indicated in iPerna fc Lazzatil l|2002ft . 
However, our analysis presented in this paper on a sample 
four times the size indicates that GRB host galaxies have 
typically larger gas-to-dust ratios than those of the SMC, 
LMC and MW. This is partly as a result of the larger sam- 
ple in this paper, but is also due to the inclusion of the UV 
and red data in the SED, which has improved the determi- 
nation of the rest frame extinction. In two cases our choices 
of SED e poch in this paper ha ve also improved the analysis, 
where in ISchadv et al.l (|2007t ) the epochs of our SEDs for 
GRB 050525A and GRB 050802 were at times when Swift 
data were poorly sampled. Nevertheless, the host galaxy 
Nh.x and Av values determined in this paper are consis- 
tent at 9 0% confidence with t he best-fit host Nh,x and Av 
values in ISchadv et all (|2007T ). 

In the rest of this section we shall explore the reasons 
that could account for the relatively large gas-to-dust ratios 
measured in GRB host galaxies. In section [6TT1 we investigate 
the effect that the GRB has on its surrounding environment, 
in section [6721 we look at the differences in the regions of gas 
probed by X-ray and optical observations, and in section [6731 



we explore the effect that the metallicity of GRB host galax- 
ies has on the Nh,x/Av ratios. 

6.1 Effect of GRB on Local Environment 

To investigate how the GRB prompt emission alters its sur- 
rounding environment, and in particular how it affects X-ray 
and optical observations of the afterglow, Perna & Lazzati 
(2002; here onwards PL02) simulated the photo-ionisation 
and dust destruction caused by a GRB within a molecular 
cloud. The more dust in the line-of-sight that is destroyed, 
and the more gas that is photo-ionised, the smaller the mea- 
sured values of Av and Nh,x, respectively. Differences in 
the efficiency of the dust destruction and photo-ionisation 
can bring about an overall change in the Nh,x/Av ratio 
measured before and after the GRB event. 

PL02 modelled the effect that a GRB would have 
on its surrounding environment when embedded within a 
Galactic-like molecular cloud with a column density of Nh = 

10 22 cm -2 and Av= 4.5 mag, and varied the particle den- 
sity by changing the radius of the cloud from R = 10 18 cm 
to R = 10 20 cm. They then simulated the soft X-ray absorp- 
tion and visual extinction that would be measured with time 
along the line-of-sight to the GRB during its prompt emis- 
sion phase, as the high energy radiation photo-ionised the 
gas and destroyed the dust within the molecular cloud. They 
found that in a large and diffuse region, photo-ionisation is 
more efficient than dust destruction, and as a result, the 
Nh,x/Av ratio measured after the GRB would be smaller 
than the value prior to the GRB. On the other hand, as the 
density becomes larger and the region more compact, dust 
destruction gradually becomes more efficient with respect 
to the photo-ionisation, and the Nh,x/Av ratio measured 
after the GRB is thus larger than the initial value. They re- 
peated their work using a GRB with a softer spectrum, and 
found a similar result but with the cross-over, when dust de- 
struction becomes more efficient than the photo-ionisation, 
occurring at a lower circumburst density. PL02 also investi- 
gated how an increase by a factor of 50 in the density within 
a cloud of radius R — 10 19 cm would effect the X-ray ab- 
sorption and optical extinction, and they found that there 
was a greater difference in the efficiency in the dust destruc- 
tion and photo-ionisation processes, where the former was 
the more effective. 

The relatively large gas-to-dust ratios measured in 
our sample of GRBs could, therefore, be a conse- 
quence of them being embedded in dense molecular 
clouds with column densities on the order of Nh = 

10 23 cm" 2 , which has already been suggested by several au- 



thors (e.g. iGalama fc Wiiersl l200ll; iReichart fc Pried [20021 : 



ICampana et al.ll2006bl ; IVergani et al.ll2004D . If dust destruc- 
tion is the cause for the large Nh,x/Av ratios, then the 
range observed in Av and Nh,x values for GRB host galax- 
ies could result from differences in the initial column densi- 
ties and sizes of the molecular cloud. A further parameter to 
consider is the location of the GRB within the cloud, where 
GRBs located closer to the outer edges nearest to the ob- 
server would destroy and photo-ionise a greater fraction of 
dust and gas along the line-of-sight than a GRB embedded 
deeper within the molecular cloud. 

One consequence of dust destruction by the GRB is the 
colour evolution of the afterglow as the dust is destroyed. 
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This is because smaller dust grains are destroyed more eas- 
ily than larger ones, such that the opacity to blue light will 
decrease sooner than the opacity to red light. However, in 
order to observe the spectral evolution caused by dust de- 
struction, multi-wavelength observations in the optical and 
infrared range are required during the dust destruction pe- 
riod. During the first ~ 600 s of the UVOT observing se- 
quence followed when a GRB occurs, only two filters are 
used. The multi-wavelength observations that are necessary 
to observe any spectral evolution in the UV-NIR afterglow 
therefore only begin ~ 600 s after the BAT trigger, by 
which time the high energy emission resp onsible for the dust 
destruction is over (jFruchter et aljfeoOll ). Robotic ground- 
based telescopes , such as the Rapid Eye Mount (REM; 
IZerbi et al1l2004h . and the RAPid Telescope fo r Optical Re- 
sponse (RAPTOR; ISalgado fc McDavidl200sh . can be tak- 
ing multi- wavelength data of a GRB a me re ~ 20 s after 
the GRB trigger (e.g. lAntonelli et al.l[2008h . However, thus 
far there have not been any clear instances of GRB light 
curves with observed colour evolution resulting from dust 
destruction. It may be possible to detect evidence of dust 
destruction in single filter observations by the increase in the 
observed flux within that filter that should observed as the 
dust opacity decreases. However, in a detailed analysis on 
the optical early-rise of a sample of six GRBs, five of which 
are i n the GRB sample studied in this paper, lOates et al.l 
(e.g. 120091 ) found no evidence of the brightening of the af- 
terglow to be the result of dust destruction for any of the 
GRBs in their sample. 

It is likely that even at 20 s after the prompt emis- 
sion, the bulk of the dust destruction has already occurred 
(PL02) , making such direct observations of dust destruction 
highly challenging for the current generation of telescopes 
and satellites. Therefore, another way of identifying a re- 
cent period of dust destruction is in the grain size distribu- 
tion in the local environment of the GRB. A consequence of 
dust destruction is a grain size distribution skewed towards 
larger grains due to the preferential destruction of small dust 
grains, which will result in a flattening of the extinction law 
llPerna et al. I l2003t ). Although large dust grains may be bro- 
ken down into smaller grains, these would, subsequently, be 
shattered promptly by the G RB, keeping the p opulation of 
small dust grains small (e.g. iPerna et al.ll2003h . The best- 
fit extinction laws to our sample of GRBs should therefore 
provide an indication of the grain size distribution of the 
extinguishing dust within the GRB environment. The fact 
that 56% of GRBs with a well-constrained host galaxy ex- 
tinction law are best-fit by the SMC extinction law (see sec- 
tion |4]2]), which has the steepest UV extinction of all the 
extinction curves fitted to our data, indicates that there is 
an abundance of small dust grains in the GRBs' surround- 
ing environments. The flattest extinction law that we fitted 
to our data was the MW extinciton law, which also has the 
most prominent 2175A absorption feature. It is, therefore, 
possible that the typically poorer fits provided by the MW 
model are due to the absence of the 21 75A absorption fea- 
ture in GRB SEDs, rather than a poor agreement between 
the slope of the MW extinction law and the GRB hosts' ex- 
tinction law. However, the LMC extinction law is natter than 
the SMC extinction law, and has a relatively weak 2I75A 
feature, and yet it was the best-fit model to only 22% of the 
sub-sample of GRBs with well-constrained power-law spec- 



tral fits. The fact that more than two times as many GRB 
SEDs are better fit with an SMC extinction law than with 
an LMC extinction law for the host galaxy suggests that 
there remains an abundance of small, UV absorbing dust 
grains in the surrounding environments of GRBs. The dust 
responsible for the UV and optical extinction must, there- 
fore, lie in regions of the GRB host galaxy that have not 
been subjected to significant amounts of dust destruction. 

6.2 Comparison between Nhi/Av and Nh,x/Av 

Nh,x provides a measurement of the optical depth of metals 
in the line-of-sight to the GRB, and in comparing the neu- 
tral hydroge n column density, Nh i, with Nh,x for a sample 
of 17 GRBs, IWatson et all (|2007h found Nhi to be system- 
atically smaller than Nh.x ■ They ascribed this to X-ray ab- 
sorption and Lya absorption observations probing different 
regions of gas. Estimates of the distance of the UV/optical 
absorbing dust from GRBs range from 100 pc to 1.7 kpc 
IjProchaska et al . 20061: IVreeswiik et al.ll2007l ). whereas there 
are likely to be partially ionised medium weight metals ab- 
sorbin g the soft X-ray emis sion within a few parsecs of the 
GRB (|Fruchter et al.ll200lT ). Although this does not neces- 
sarily mean that all the neutral hydrogen out to 100 pc is 
ionised, it does suggest that all the neutral hydrogen within 
the molecular cloud surrounding the GRB has been ionised. 
So whereas Nhi will probe the host galaxy ISM outside 
of the molecular cloud, Nh,x will, in addition, probe the 
gas within the molecular cloud, where densities of partially 
ionised oxygen and other medium weight metals remain rel- 
atively high. These differences in the gas probed by Nhi and 
Nh.x measurements provide information on the conditions 
of the environment at varying radii from the GRB. Having 
compared Ay to Nh,x, we now look at the relation between 
Ay and Nhi to investigate the relative location of the re- 
gions of dust and gas probed by these two measurements. 

We took Nhi values reported in the literature for all 
GRBs that overlapped with our sample (Jakobsson et al. 
2006 and references therein), of which there are eight (see 
Table EJ. Fig. H] shows Ay from the SMC (top panel), LMC 
(middle panel) and MW (bottom panel) spectral model fits 
plotted against Nhi, in log — log space. The dashed line 
in each panel, from top to bottom, repr esent th e Nhi I Ay 
ratios for the SMC dMartin et aflll989l), LMC (iKoornnee: 
ll982tlFitzpatricldll985h and MW (Predehl fc Schmittlll9<T ~ 
respectively. 

For each GRB plotted in Fig. |9j Nhi is typically an 
order of magnitude smaller than Nh,x as was noted in 
IWatson et all (120071 ). and the GRB N H i/A v ratios span at 
least an order of magnitude to each side of the Milky Way 
and Magellanic Cloud Nhi/ Ay ratios. This is in contrast to 
the GRB Nh,x/Av ratios, which were mostly larger than 
those of the Milky Way and Magellanic Clouds. A spearman 
rank test between Nhi and the best- fit Ay from the SMC, 
the LMC and MW models gives the spearman coefficients 
0.34, 0.25 and 0.09 respectively, indicating that there is only 
a weak correlation, if any, between these two parameters. 

Where the GRB data points lie relative to the dashed 
lines in Fig. [9] is the result of two competing effects. Host 
galaxy metallicities that are smaller than the Magellanic 
Clouds, as well as any significant dust destruction caused 
by the GRB, will move the data points downwards, below 
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Table 3. Subsamplc of GRBs with Nhi and metallicity, [M/H], 
measurements available in the literature. 



GRB 


log N HI 


[M/H] 


050319 


20.9 iO.2 1 




050730 


22.1 iO.l 1 


-2.26±0.14 3 


050820A 


21.1 iO.l 1 


-0.63±0.11 3 


050922C 


21.6±0.1 1 


-2.03±0.14 3 


060124 


18.5 ±0.5 2 




060206 


20.9 iO.l 1 


-0.85±0.18 3 


060418 




-1.65±1.00 3 


060526 


20.0 iO.2 1 


-1.09±0.24 4 


060607A 


< 19. 5 1 




060714 


21.80 ±0.10 2 




070110 


21.70 ±0.10 2 




070411 


19.30 ±0.30 2 





1 IJakobsson et al.l j2006ah 



Fvnbo et al.l |2009|) 



3 iProchaska et al 
4 iThoene et al.l 



,2007) 
[ 20081) 



the dashed lines, and counteracting this effect will be the 
amount of photo-ionisation of hydrogen in the surrounding 
environment of the GRB that will move the data towards 
the left of Fig. [9] The mean logarithmic metallicity of those 
GRBs plotted in Fig[H]with known metallicity is {[M/H]) = 
— 1.37 (i.e. 0.04 Zq), which is almost 1.0 dex smaller than 
the SMC, in which case we would expect the majority of 
the data points in Fig. [5]to lie below the dashed lines. The 
fact that the data points are fairly evenly distributed about 
the dashed lines, therefore, indicates that a greater volume 
of gas has been photo-ionised by the GRB, thus reducing 
Nhi, than the volume of dust destroyed, which reduces 
Ay. This is consistent with the analysis of high-resolution 
spectroscopic data, which indicates that GRBs can photo- 
ionise all gas within the molecular cloud surrou nding the 
GRB l|Prochaska et al.ll2007l ; IVreeswiik et al.ll2007f ). whereas 
a GRB is only expected to fully destr oy dust out to a few 
parsecs at most ijPerna fc Lazz ati 2002). The t ypically larger 
values of Nh.x compared to Nhi discussed bv lWatson et al.l 
l|2007l ). and the distribution in Nhi / 'Av plotted in Fig. [9] 
therefore, suggest that measurements of Nh,x and Av probe 
regions of gas and dust within the molecular cloud, much 
closer to the GRB than measurements of Nhi- 



6.3 GRB host galaxies and Dwarf Irregulars 

It has been noted that for a large range of galaxy types, 
the gas-to-dust ratio of galaxies is inversely proportional 
to the metallicity, down to the most metal-poor systems 
ijDraine et al.ll2007l ). The range in column density to ex- 
tinction ratio in GRB host galaxies may, therefore, be ac- 
counted for if those host galaxies with lower metallicities 
have larger Nh,x/Av ratios. Evidence of such a correla- 
tion in our GRB sample is indicated in Fig. 1101 Here we 
have plotted Nh,x/Av against the host metallicity, [M/H], 
for a subsample of GRBs (solid circles) that have an esti- 
mate of the host metallicity as well as a host galaxy soft 
X-ray absorption system and/or a dust extinction system 
detected with 90% confidence. The top, middle and bot- 
tom panels show the Nh,x/Av ratios determined from the 
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Figure 9. Logarithmic host galaxy Ay against logarithmic Nhi 
for a subsample of GRBs with Nhi mea surements availab l e in the 
literature. Nhi values are taken from IJakobsson et al" 
and references therein. The Ay values are the results from our 
spectral analysis from the SMC model (top panel), the LMC 
model (middle panel) and the MW model (bottom panel). The 
dashed curves are the Nhi/ Ay ratios for each corresponding en- 
vironment. 



SMC, LMC and MW spectral model fits, respectively. Those 
GRBs in our sample with host metallicity measurements 
are listed in Table [3] along with their metallicity and corre- 
sponding reference. Although listed in Table GRB 060206 
and GRB 060526 are not included in Fig. 1101 since they only 
have upper limits for both Nh,x and Ay. The open circles 
in Fig. [TD] correspond to the SMC, LMC and Milky Way 
from left to right, respectively. 

If very low metallicities are the over-riding reason 
that GRB host galaxies typically have larger Nn,x/Ay 
ratios than the Magellanic Clouds and Milky Way, then 
there should be a correlation between the metallicity and 
Nh,x I Ay ratio for the combined sample of GRB host galax- 
ies and other, more metal-rich galaxies. A spearman rank 
test between the Nh,x/ Ay ratio and the metallicity, [M/H], 
for the four GRBs shown in Fig. [10] together with the SMC, 
LMC and MW data points gives a coefficient of -0.89 with 
90% confidence for each of the spectral models. This in- 
dicates a strong anti-correlation with a high level of signifi- 
cance. The dashed line in each panel is the best-fit power-law 
to the data. 

It such a correlation is confirmed in the future with 
a greater sample of GRB host galaxies with a measured 
metallicity and Nh,x /Ay ratios, this would imply that low- 
metallicity galaxies are less efficient at forming dust from 
their metals than high-metallicity galaxies. One possible 
cause of this is an increase in supernovae dust destruction ef- 
ficiencies in low metallicity environm ents, resulting from in - 
termittent periods of star formation ( Hirashita et al.ll2002h . 
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Figure 10. Host galaxy Nh,x I Ay against metallicity, [M/H], 
for a subsamplc of five GRBs (solid circles) with [M/H] values 
available from the literature (see Table [3j and a soft X-ray ab- 
sorption system and/or dust extinction system detected with 90% 
confidence. The Nh X and Ay values are the best-fit parameters 
from the SMC (top), LMC (middle), and MW (bottom) spectral 
fits. Open circles correspond to the SMC, LMC and Milky Way 
from left to right, respectively. The dashed line is the line of best- 
fit to the SMC, LMC, Milky Way, and the four GRBs with a soft 
X-ray absorption and dust extinction system detected with 90% 
confidence. 



7 SUMMARY & CONCLUSIONS 

In this paper we have presented the results from the spectral 
analysis of 28 GRB SEDs. We measured the equivalent neu- 
tral hydrogen column density and visual extinction at the 
host galaxy, and found 79% of the GRBs in our sample to 
have a detectable soft X-ray absorption system in the host 
galaxy, and 71% to have a detectable visual dust-extinction 
system. Using the measured Nh,x/Av ratios as an indicator 
of the host galaxy gas-to-dust ratio, we find that GRB host 
galaxies have gas-to-dust ratios that are typically larger than 
those measured in the Milky Way and Magellanic Clouds by 
up to two orders of magnitude. We have investigated several 
possibilities that could account for the relatively large gas- 
to-dust ratios in GRB host galaxies. 

There is no evidence to suggest that the large host 
galaxy Nh.x /Ay ratios measured in our GRB sample is the 
result of any systematic error in the way that we measure 
Ay. One possibility is that dust destruction by the GRB has 
reduced the visual extinction, Ay, relative to the equivalent 
neutral hydrogen column density, Nh,x- However, there are 
currently no observations that clearly show the early time 
colour evolution expected from dust destruction. Although 
such observations are limited by the quality and prompt- 
ness of the data, we also found that the majority of our 
sample had host dust properties best-fit by the UV steep, 
SMC extinction law, indicating an abundance of small dust 



grains in the GRB surrounding environment. In the event 
of a significant phase of dust destruction, a grey extinction 
law should be observed, where the differential change in ex- 
tinction from UV to NIR energy bands is small. The dust 
probed by our Ay measurements must, therefore, lie in re- 
gions of the GRB host galaxy that have not been subjected 
to significant amounts of dust destruction. 

For a subset of eight GRBs we were also able to study 
how the neutral hydrogen column density, Nhi, compared 
with Ay, and we found Nhi/ Ay to extend to both larger 
and smaller values than those of the Magellanic Clouds and 
the Milky Way by up to an order of magnitude. The distribu- 
tion in Nhi/ Ay can be accounted for by the competing ef- 
fects that alter the values of Nhi and Ay. Firstly, differences 
in the host galaxy metallicities and in the amount of dust de- 
stroyed by the GRB will affect the value of Ay. On the other 
hand, the value of Nhi will be dependent on the amount of 
photo-ionised hydrogen along the line-of-sight to the GRB. 
The mean logarithmic metallicity of the GRB sample with 
both Nhi and Ay measurements is almost 1.0 dex smaller 
than that of the SMC (0.04 Zq), and we would therefore 
expect the GRB host galaxy Nhi/ Ay ratio to be signifi- 
cantly smaller than the SMC Nhi/ Ay ratio. The roughly 
even number of GRBs with smaller and larger Nhi /Ay ra- 
tios than the Magellanic Clouds and Milky Way therefore 
implies that the level of photo-ionised hydrogen along the 
line-of-sight to the GRB is greater than the fraction of dust 
destroyed by the GRB. This would suggest that measure- 
ments of Nh,x and Ay probe regions of dust and gas much 
closer to the GRB than Nhi- 

It has been suggested that differences in the gas-to-dust 
ratios in galaxies of different types are correlated with the 
metallicity of the galaxy (e.g. iDraine et al.ll2007r i. whereby 
smaller metallicity systems have larger gas-to-dust ratios. 
From a subsample of four GRBs with measured metallic- 
ity and a soft X-ray absorption and visual extinction sys- 
tem detected with 90% confidence, together with the Small 
and Large Magellanic Clouds and Milky Way, we found a 
strong negative correlation between the Nh,x/ Ay ratio and 
the metallicity, [M/H]. The spearman rank coefficient was 
-0.89 with 90% confidence. The large Nh.x/ Ay ratios mea- 
sured in GRB host galaxies could, therefore, be an indication 
of their very low, although broad, range of metallicities. A 
greater sample of GRB hosts with measured metallicities 
are needed to verify such a correlation, which if confirmed 
would suggest that low-metallicity environments are less effi- 
cient at forming dust from their metals than high-metallicity 
galaxies. 
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Table 4. Results from simultaneous UV/optical and X-ray spectral fits for the SMC, LMC and 
MW dust-extinction law models, for both a power-law (pow) and a broken power-law (bknp) 
continuum. The third and fourth columns give the host galaxy equivalent column density and 
visual extinction, the fifth column gives the break energy for the broken power-law spectral models, 
the x 2 an d degree of freedom (dof) of the fit are given in the sixth column, and the seventh gives 
the null hypothesis probability. 



GRB Model N H}X A v E bk x 2 (dof) Null Hypothesis 

10 21 cm" 2 (mag) (Hz) Probability 



050318 SMC/pow l-40Io 40 °- 53 -o.06 ' 108 ( 88 ) °- 068 

LMC/pow l-9ltg:t? a78 ±oio9 ' 113 ( 8S ) °- 039 

MW/pow 1.67j:o_4g 0.881°;" - 147 (88) 8.3c-05 

SMC/bknp 1.4lj;°;|| 0.54±°;^ 8.616 108 (87) 0.059 

LMC/bknp 1.92lg;|| 0.79j;°;^ 8.971 113 (87) 0.034 

MW/bknp 1.67±g;|| 0.88t°;J 2 9.158 147 (87) 6.2c-05 



050319 SMC/pow < 3.40 0.07Iq £3 - 81 (80) 0.448 

LMC/pow < 4.36 0.12+°'°;! 80 (80) 0.492 

MW/pow < 5.25 0.21±g;gf - 77 (80) 0.588 

SMC/bknp < 4.16 < 0.09 2.541 71 (79) 0.714 

LMC/bknp < 4.45 < 0.21 2.566 71 (79) 0.722 

MW/bknp < 5.01 < 0.35 2.636 70 (79) 0.755 



050525A SMC/pow 2.25I l ; 4 g 0.06 ± 0.01 50 (37) 0.069 

LMC/pow 2.24±° ) ; 4 ( 1 0.07 ± 0.02 - 56 (37) 0.024 

MW/pow 2.15±g;|J 0.06l°;° 2 - 68 (37) 0.001 

SMC/bknp 2.96 + .° ) ; 4 g 0.16t°;o2 °- 017 46 ( 36 ) °- 119 

LMC/bknp 3.12 ±0.53 0.21 ± 0.03 0.022 57 (36) 0.014 

MW/bknp 3.82 ±0.73 0.10 ±0.04 0.724 84 (36) 1.0c-05 



050730 SMC/pow 14.88+. 2 ; 2 3 al6 -o'.o1 163 ( 133 ) °- 040 

LMC/pow 15.51±|;^g 0.22 ±0.03 163 (133) 0.039 

MW/pow 16.20+.2 0.30 ± 0.05 164 (133) 0.036 

SMC/bknp 17.79l 2 ;^ °- 23 tao3 °- 001 159 ( 132 ) °- 055 

LMC/bknp 18.31^2.01 a31 ±o:o5 a001 159 ( 132 ) °- 053 

MW/bknp 18.14t 2 ;^ 0.39^°;^ 0.001 160 (132) 0.047 



050802 SMC/pow 1.28±l\H 0.06 ± 0.02 - 90 (69) 0.043 

LMC/pow 1-43^0 56 0.10 ± 0.03 89 (69) 0.056 

MW/pow 1.74+.°;|* 0.19 ±0.06 84 (69) 0.101 

SMC/bknp 1.56±q; ™ 8 0.05+°;q 2 3.150 81 (68) 0.138 

LMC/bknp 1.67to 5° 0.08+. ,; , 2 3 - 168 80 ( 68 ) °- 155 

MW/bknp 1.89±g;|f 0.15+°;°^ 3.233 77 (68) 0.212 



050820A SMC/pow < 0.46 0.181°;°} 224 (139) 6.2e-06 

LMC/pow < 0.67 0.29^°;°^ 208 (139) 1.3c-04 

MW/pow < 1.45 0.43t°;°4 193 (139) 0.002 

SMC/bknp 5.07^0,11 0.14 ±0.03 0.209 144 (138) 0.350 

LMC/bknp 5.04±{™ 0.23 ± 0.04 0.143 142 (138) 0.386 

MW/bknp 5.35±i"_2i 0- 32 ±o.o6 0139 143 ( 138 ) °- 365 



050922C SMC/pow 1.62j;g-|| 0.07 ± 0.02 - 36 (48) 0.902 

LMC/pow l-73+°;^° 0.11+°;o3 3 5 (48) 0.915 

MW/pow l-85l°;° 2 0.16to°l 35 (48) 0.925 

SMC/bknp 3.03^1% 0.14+. ,; 02 0.005 34 (47) 0.915 

LMC/bknp 3.13±^;|? 0.21+°;°^ 0.005 34 (47) 0.926 

MW/bknp 3.581};°* 0-^to° Q e 5 0.007 37 (47) 0.851 
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Table 4. Continued 



GRB Model N H>X A v E bk x 2 (dof) Null Hypothesis 

10 21 cm" 2 (mag) (Hz) Probability 



051109A SMC/pow < 3.64 < 0.14 239 (145) 1.3e-06 

LMC/pow < 3.64 < 0.19 - 239 (145) 1.3e-06 

MW/pow < 3.64 < 0.25 - 239 (145) 1.3e-06 

SMC/bknp 8.7lj;^ < 0.10 0.090 162 (144) 0.141 

LMC/bknp 8.75±i;|| < 0.13 0.092 162 (144) 0.141 

MW/bknp 8.80 ±1.55 < 0.17 0.096 162 (144) 0.140 



060124 SMC/pow 2.60I° |j 0.08 ± 0.03 182 (115) 6.6e-05 

LMC/pow 3.24tg'f 6 0.17 ±0.04 174 (115) 3.0c-04 

MW/pow 5.83±J;J| 0.56±° °| - 144 (115) 0.034 

SMC/bknp 9Al±{-* 7 7 0.13±%%l 0.233 116 (114) 0.439 

LMC/bknp 9.39+^33 0.22 ± 0.06 0.156 114 (114) 0.483 

MW/bknp 9-68iJ;3i 0.52^°;^ 0.060 113 (114) 0.498 



060206 SMC/pow < 6.72 < 0.04 - 97 (53) 2.0e-04 

LMC/pow < 6.72 < 0.05 97 (53) 2.0c-04 

MW/pow < 6.72 < 0.05 97 (53) 2.0c-04 

SMC/bknp 13.65±2'.?g < °- 18 °- 618 60 ( 52 ) °- 218 

LMC/bknp 14. 06^3.92 < °- 05 °- 604 60 ( 52 ) °- 218 

MW/bknp 13.69^2.83 < °- 22 °- 617 60 ( 52 ) °- 218 



060418 SMC/pow 3.60ti ll < om 41 ( 24 ) °- 017 

LMC/pow 3.53+^;| 2 < 0.09 - 42 (24) 0.014 

MW/pow 3.44tj j° < 0.10 42 (24) 0.013 

SMC/bknp 4.02± 2 ;3 2 0.09t°;o 2 0.002 21 (23) 0.591 

LMC/bknp 4.24+ 2 ; 27 i 0.12+°;° 2 0.001 23 (23) 0.461 

MW/bknp 3.65±5:li °- 08 ±ao2 °- 001 29 ( 23 ) °- 186 



060502A SMC/pow 3.42t° j° 0.5li°;^ 33 (31) 0.368 

LMC/pow 4.03±g;|| 0.74± \l 35 (31) 0.299 

MW/pow 3.88±J;^| 0.79±g;if 55 (31) 0.005 

SMC/bknp 5.07lJ;|j 0.50±°;Jg 0.031 30 (30) 0.467 

LMC/bknp 4.10+Q °g °- 76 lao6 °- 002 34 ( 30 ) °- 297 

MW/bknp 6.18lJ;^ 0.50tg:?g 0.432 57 (30) 0.002 



060512 SMC/pow < 0.85 0.47 ± 0.05 84 (23) 7.4e-09 

LMC/pow < 0.85 0.56 ±0.06 - 99 (23) 2.1c-ll 

MW/pow < 0.82 0.67 ±0.08 - 122 (23) 2.0c-15 

SMC/bknp < 1.74 0.66 ± 0.09 0.007 77 (22) 4.8e-08 

LMC/bknp < 1.79 0.79 ±0.13 0.011 99 (22) 1.0c-ll 

MW/bknp < 2.03 0.95 ±0.15 0.015 131 (22) 2.1c-17 



060526 SMC/pow < 45.39 < 0.07 9 (8) 0.344 

LMC/pow < 45.61 < 0.10 - 9 (8) 0.344 

MW/pow < 45.67 < 0.21 9 (8) 0.345 

SMC/bknp < 47.14 < 0.16 0.002 7 (7) 0.476 

LMC/bknp 13.23±l 2 gl 5 0.10 ±0.04 0.002 6 (7) 0.486 

MW/bknp 14.691^4^ °- 18 ±ao9 °- 002 6 ( 7 ) °- 554 
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Table 4. Continued 



GRB Model N H>X Ay E bk x 2 (dof) Null Hypothesis 

10 21 cm" 2 (mag) (Hz) Probability 



060605 SMC/pow 7.80+ 2 '™ 0-251°;°° 62 (68) 0.670 

LMC/pow 7.79+^;^ 0.32l°;o7 62 (68) 0.673 

MW/pow 6.47±| || 0.35±^;g| 63 (68) 0.662 

SMC/bknp 7.70l 2 ;™ 0.24+.°;°® 4.892 62 (67) 0.662 

LMC/bknp 7.61+J'g 2 0.31+ ( ( 5;°g 4.971 62 (67) 0.664 

MW/bknp 6.13±2'44 0.3ll°;Jg 3.947 62 (67) 0.666 



060607A SMC/pow 6.19+^;™ < 0.15 - 116 (117) 0.506 

LMC/pow 6.16lJ ®3 < 0.21 - 116 (117) 0.499 

MW/pow 5.961^2 < 0.26 - 117 (117) 0.489 

SMC/bknp 6.27±^'|| < 0.11 9.787 116 (116) 0.479 

LMC/bknp 6.22+^ < 0.14 9.786 116 (116) 0.472 

MW/bknp ^ < 0.14 9.832 117 (116) 0.462 



060714 SMC/pow 5.98+3;°2 0.46+.°;^ 20 (17) 0.284 

LMC/pow 6.48^91 0-64l°; 2 ^ - 20 (17) 0.268 

MW/pow < 20.04 0.79t°;|° - 22 (17) 0.196 

SMC/bknp < 17.07 < 0.91 2.341 19 (16) 0.283 

LMC/bknp < 17.86 < 1.28 2.241 19 (16) 0.261 

MW/bknp < 16.87 < 1.66 2.121 20 (16) 0.211 



060729 SMC/pow 0.80 ±0.09 0.03 ± 0.01 - 184 (178) 0.367 

LMC/pow 0.82 ± 0.09 0.04 ± 0.02 - 183 (178) 0.373 

MW/pow 0.83l° J° 0.06 ± 0.03 184 (178) 0.354 

SMC/bknp 1-09+.°;" 0.13+°;^ 0.006 176 (177) 0.496 

LMC/bknp LlOl^Jg 0.18l°;o2 °- 006 17 8 (177) 0.465 

MW/bknp 0.971°;°?, °- 15 to° l 0.004 179 (177) 0.453 



060904B SMC/pow IM+q^I 0.06 ± 0.02 - 83 (46) 7.3c-04 

LMC/pow 1.85to 37 0.08 ± 0.03 83 (46) 6.3c-04 

MW/pow 1-84+°;^ 0.09 ± 0.04 85 (46) 3.7e-04 

SMC/bknp 3.72±°;^ 0.12+°;°^ 0.174 74 (45) 0.004 

LMC/bknp 3.65+_°,;£g 0.17 ±0.06 0.129 75 (45) 0.003 

MW/bknp 4.05+°;^ 0.14+°;°* 0.392 80 (45) 0.001 



060908 SMC/pow < 8.74 < 0.17 - 15 (8) 0.061 

LMC/pow < 8.95 < 0.38 15 (8) 0.060 

MW/pow < 8.71 < 0.22 - 15 (8) 0.061 

SMC/bknp < 25.95 < 0.21 1.287 12 (7) 0.111 

LMC/bknp < 13.77 < 0.26 1.212 12 (7) 0.109 

MW/bknp < 25.82 < 0.18 1.288 12 (7) 0.111 



060912 SMC/pow 3.23i°-g2 O- 44 !"'.^ - 60 ( 42 ) O- 037 

LMC/pow 3.52t°'|^ 0.621°;^ 59 (42) 0.040 

MW/pow 0.8lt°;^ 61 (42) 0.027 

SMC/bknp 3.23+°j° 0.44±°;i° 5.975 60 (41) 0.029 

LMC/bknp 3.53l°;| 2 0.62+°;°® 8.717 59 (41) 0.031 

MW/bknp 3.70+.°;!° °- 82 t° °l 7.695 61 (41) 0.021 
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Table 4. Continued 



GRB Model N H>X A v E bk x 2 (dof) Null Hypothesis 

10 21 cm" 2 (mag) (Hz) Probability 



061007 SMC/pow 0.45 ± 0.01 380 (274) 2.2e-05 

LMC/pow 5.44+°; 2 J 0.75 ±0.02 325 (274) 0.019 

MW/pow 6.82 ± 0.23 1.26 ± 0.03 - 1157 (274) 0.0c+00 

SMC/bknp 5.10+°H 0.53±%%\ 0.004 337 (273) 0.005 

LMC/bknp 5.74 ± 0.21 0.82+j^J 0.003 315 (273) 0.040 

MW/bknp 6.79 ±0.26 1.26 ±0.03 9.992 1157 (273) 0.0c+00 



061121 SMC/pow 3.64I°;~ 0.34 ± 0.03 - 176 (131) 0.005 

LMC/pow 4.04+°;^ 0.49 ±0.04 167 (131) 0.018 

MW/pow 4A6+°'l° 0.71 ± 0.06 163 (131) 0.031 

SMC/bknp 3.87+.°;^ 0.28 ± 0.03 2.738 139 (130) 0.272 

LMC/bknp 4.13l°;le 0.40 ± 0.04 2.829 138 (130) 0.305 

MW/bknp 4.32+.° l ;g 0.55 ± 0.07 2.865 139 (130) 0.270 



061126 SMC/pow 3.14+° 3 2 0.13±[^ - 190 (133) 8.9e-04 

LMC/pow 3.35+°;*{ 0.21+°;°* 186 (133) 0.002 

MW/pow 3.50±°,'43 0.28lg;°g 186 (133) 0.002 

SMC/bknp 5.95±° ™ 0.10 ±0.04 0.135 144 (132) 0.225 

LMC/bknp 5.97±° l* 0.14t°;°^ 0.122 145 (132) 0.209 

MW/bknp 0.ut ° 7 7 0.170 149 (132) 0.153 



070110 SMC/pow < 3.29 0.29+.°,;°° 55 (50) 0.306 

LMC/pow < 4.23 0.44t° ™ 51 (50) 0.432 

MW/pow < 5.34 0.64^^2 50 (50) 0.472 

SMC/bknp < 3.48 0.23±° ° 6 5 2.408 45 (49) 0.656 

LMC/bknp < 4.07 0.341°;°^ 2.530 44 (49) 0.692 

MW/bknp < 4.69 0.49±[J;ig 2.828 44 (49) 0.662 



070318 SMC/pow 8.52t°;*3 0.52 ± 0.02 69 (44) 0.010 

LMC/pow 8.96 ±0.77 0.73 ± 0.03 94 (44) 1.5c-05 

MW/pow 9.49 ± 0.81 1.05 ± 0.04 - 213 (44) 4.8c-24 

SMC/bknp S^to.el °- 59 ±o.o6 °- 003 65 ( 43 ) °- 017 

LMC/bknp 9.21 ±1.14 0.82 ± 0.05 0.003 86 (43) l.lc-04 

MW/bknp 9.21 ±1.13 0.82 ± 0.05 0.003 86 (43) l.lc-04 



070411 SMC/pow < 32.05 < 0.21 - 24 (17) 0.117 

LMC/pow < 31.93 < 0.30 24 (17) 0.117 

MW/pow < 31.93 < 0.47 - 24 (17) 0.117 

SMC/bknp < 32.81 < 0.20 5.026 24 (16) 0.089 

LMC/bknp < 32.68 < 0.28 5.026 24 (16) 0.089 

MW/bknp < 32.71 < 0.44 5.026 24 (16) 0.089 



070529 SMC/pow 14.12I°;|2 < 0.52 20 (25) 0.728 

LMC/pow 14.23^;^ < 0.68 20 (25) 0.729 

MW/pow 13.76±5.?9 < 0.84 20 (25) 0.724 

SMC/bknp 13.52+^-j;? < 0.45 3.734 18 (24) 0.780 

LMC/bknp 13.61±f*l < 0.58 3.698 18 (24) 0.780 

MW/bknp 13.41^; 2 2 < 0.70 3.713 18 (24) 0.781 



